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ABSTRACT 


The contribution of pulsars to the y-ray flux from the galactic plane 
is examined using data from the most recent pulsar surveys. It is assumed 
that pulsar y-rays are produced by curvature radiation from relativistic 
particles above the polar cap and attemiated by pair production in the 
strong magnetic and electric fields. Assuming that all pulsars produce 
y-rays in this way, their luminosities can be predicted as a function of 
period and magnetic field strength. The distribution of pulsars in the 

Galaxy is determined from data on 328 pulsars detected in three surveys. 

The Z and R distributions are very sensitive to the mean electron density 
in the galactic plane, as are the total number and birthrate of pulsars in 
the Galaxy. 

The local y-ray production spectrum for pulsars is steep above 100 MeV 
and is similar to the bremsstrahlung and Compton spectra. Longitude profiles 
of pulsar y-ray flux are calculated for different values of the mean electron 
density. Because of the large center to anticenter flux ratio, pulsars 
contribute twice as much to the total flux toward the galactic center as they 
do toward the anticenter. The latitude profile is narrow, due to the fact 
that short period pulsars, which have the highest y-ray luminosities, also 
have the smallest scale heights. The largest sources of uncertainty in the 
size of the pulsar contribution are the value of the mean interstellar electron 
density, the turnover in the pulsar radio luminosity function, and the 
average pulsar magnetic field strength. A present estimate is 
that pulsars contribute from 15-20 percent of the total flux of y-rays from 
the galactic plane. 


I. INTRODUCTION 


The plane of our galaxy Is the dominant source of high energy (> 100 MeV) 
Y-rays detected with the SAS-2 (Fichtel et al. 1975, Hartman et al, 1979) and 
COS-B (Paul et al. 1978, Mayer-Hasselwander et al, 1979) satellites. These 
observations have revealed, in some detail, the longitude and latitude 
dependent structure of the galactic y-radiation. The large-scale emission 
shows a pronounced peak toward the galactic center between 330° and 
4*- 50° and is relatively weak in the outer galaxy. The latitude profile 
is fairly narrow and falls off significantly within 10° of b “ 0°. On a 
finer scale, the emission shows spatial fluctuations, some of which may 
correspond to local or spiral arm features and all of which a»e smaller 
than the 1° resolution limit of the detectors. Two of these stand out strongly 
in the longitude profile and have been identified* through their timing 
signature, as the Crab and Vela pulsars. A number of other “point sources" 
of emission have been identified by COS-B (Wills et al. 1980), all of which have 
so far eluded definite identification with galactic sources at other v^avelengths, 
although models and identifications have been proposed for some of the sources. 

Diffuse emission processes which involve interactions between high 
energy cosmic rays and interstellar matter are considered to be a major 
source of galactic Y-rays. These processes include the decay of neutral 
pions produced in collisions between cosmic ray nucleons and interstellar 
gas nuclei, bremsstrahlung from cosmic ray electrons in the Coulomb fields of 
nuclei, and Compton scattering of microwave background and starlight photons 
from the cosmic ray electrons. The y-fay emission, therefore, has great 
potential for revealing information on the origin and distribution of cosmic 
rays in the Galaxy and on the distribution of matter. Models for the galactic 
emission from diffuse processes using known distributions of CO and HI 
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(stecker et al. 1975, Kniffen et al. 1977) require an increase in 
the cosmic ray density in the inner galaxy to explain the peak 
in y-ray emission toward the galactic center, even with a large increase in 
molecular hydrogen density in the inner galaxy. This strongly suggests a 
galactic origin for most cosmic ray nuclei (Stecker 1975). 

This picture, however, is complicated by the largely undetermined con- 
tribution from galactic point sources. Estimates have been made of the point 
source contribution using data on the COS-B sources (Protheroe et al. 1979, 
Bignami et al, 1978), but the estimates are uncertain due to a lack of 
knowledge of the distribution of these sources in the Galaxy. Without identi- 
fication at other wavelengths, distance determination is impossible and thus 
the luminosities of these sources are also unknown. Another difficulty with 
this type of analysis is the inability to distinguish true point sources from 
enhancements in the spatial distribution due to diffuse processes. 

At present, pulsars can give the best information on the galactic point 
source contribution. Since independent distance determinations can be made 
for pulsars via their dispersion measures, the galactic distribution of these 
objects is reasonably well known, except in the inner 3 kpc where statistics 
are poor. The Crab and Vela pulsars emit pulsed y-rays and have been well 
studied by COS-B and SAS-2 (Bennett et al. 1977, Kniffen et al. 1974). Although 
a theoretical model for the acceleration of particles to energies high enough 
to produce the observed y radiation has not been generally agreed upot}, pulsars 
are certainly good candidates for high energy y-ray sources on energetic grounds. 
Assuming that an efficient acceleration mechanism is operative, models for the 
production of y-rays by primary particle curvature radiation, including the 
attenuation by pair production, are able to account for many of the observed 
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properties of the Crab and Vela pulsars (Harding et al. 1978, Salvati and 
Massaro 1978). Theory also suggests that all pulsars should be y-ray emitters, 
with the y-ray efficiency an increasing function of age and the shortest period 
pulsars having the highest luminosities (Harding 1981, Ayasli and Ogelman 
1980). Furthermore, the predicted fluxes for the known radio pulsars 
Indicate that many are right below the sensitivity threshold of the present 
y-ray detectors. 

Previous estimates have been made of the pulsar contribution to the 
y-ray emission above 100 MeV. Higdon and Lingenfelter (1976) suggested that 
the contribution from unresolved pulsars could be as high as 40 percent, 
assuming that their y-ray luminosities are directly proportional to their 
energy loss rates and that they have the same galactic distribution as the 
CO emission. Strong et al. (1977), assuming a direct proportionality 
between y-ray and radio luminosity and using the distributions determined 
by Davies et al. (1977), concluded that the contribution would be no more 
than 5-10 percent. Neither of these estimates made use of a theoretical 
model for pulsar y-ray emission and the large discrepancy in the results 
seems to indicate large uncertainties or errors in the assumptions made. 

This paper makes a mere thorough examination of the pulsar contribution 
and of the sources of uncertainty involved. We begin with a model for pulsar 
y-ray production (§11) which can predict the luminosities as a function of 
pulsar parameters. All of the present data on radio pulsars is then used 
to determine the galactic distribution, taking into account the selection 
effects of the surveys (§III). In §IV, we estimate the pulsar contribution 
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to the local y-ray einissivity and compare it to the production rates for 
various diffuse processes. Calculated longitude and latitude profiles are 
presented in §V and compared to the SAS-2 data. We conclude In 5VI by 
discussing the major sources of uncertainty and the most reasonable present 
estimate of the pulsar contribution. 



II. PULSAR Y-RAY luminosities 
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The model v^e will use for the production of pulsar v-rays is described in 
previous work (Harding et al. 1978, Harding 1981). It assumes that primary 
particles are accelerated to high energies in the magnetospheri c electric fields 
and move along curved magnetic field lines, losing energy by generating curvature 
radiation y-rays. The y-rays which will be observed are those which escape 
conversion to electron-positron pairs in the strong pulsar magnetic field. The 
model takes into account rotation effects such as aberration, the electric field 
which is induced perpendicular to the magnetic field, and the rotation of the 
dipole field pattern in the frame of the photon. The shape of the calculated 
spectra depend on the initial energy of the primary particles, and on the 
pulsar magnetic field strength and period. A model with initial particle 
energy 1.5 X 10' ’ eV and field strength 10'* gauss gives a good fit to the 
observed spectra of both the Crab and Vela pulsars. 


A general luminosity formula was derived and found to depend on the 
period P (in seconds) and the surface magnetic field strength Bfs (in 
units of 10^* gauss) as: 

L^(> 100 MeV) = 1.2 x lO^® P'"*’ Photons s“\ (1) 

which is normalized to give the correct luminosity for the Crab. This 
formula gives a value of L (> 100 MeV) for Vela which is in good 

d 

agreement with the observed value. Predicted luminosities for 
other pulsars give fluxes which, with the exception of a few, are all 
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below the upper limits given by SAS-S (Ogelman et al, 1975). The few 
which are above the SAS-2 upper limits only exceed them by a factor of 2 
or 3, which is not a serious discrepancy, considering that the errors in 
the limits are of the same magnitude. In addition, individual pulsar 
distances and magnetic field strengths can only be determined to within 
a factor of 2 or 3. The second COS-B catalog of point sources is now 
considered complete in the region 90° < / < 270° down to a flux of 1.3 
X 10“® cm-® S-’ (Swanenburg et al. 1981). Only two pulsars with predicted 
fluxes above this value lie in that longitude range and they are both 
located at high latitudes just outside the range of the search. None of 
the other pulsars have predicted fluxes exceeding the observed diffuse 
background level in their region of the sky. It seems that, at present, 
there are no serious discrepancies between what is predicted by the model 
and what is observed. Therefore, we will adopt the formula in equation (1) 
as a good estimate, on the average, of pulsar y-ray luminosities, since 
we are concerned in this paper only with the collective properties of 
pulsar y-ray emission in the Galaxy. 
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III. DISTRIBUTION OF PULSARS IN THE GALAXY 

It is known from earlier statistical analyses that pulsars arc 
distributed much like ether Population I tracers in the Gala>^, such as 
ionized hydrogeni CO, and supernova remnants. The densities rise sharply 
toward the inner galaxy and fall off outside the solar circle. Previous 
analyses have been carried out on samples of 51 (Davies et al. 1977) 

90 (Taylor and Manchester 1979), and 224 pulsars (Manchester 1979). 

At present, there are 328 pulsars which have been detected in the 
various radio surveys. The three most sensitive searches, carried out at 
Arecibo (Hulse and Taylor 1974), Holonglo (Manchester et al. 1978) and 
NRAO (Damashek et.al. 1978), have collectively detected all the pulsars 
discovered in previous less sensitive surveys plus a substantial number 
of new ones. The entire galactic plane and all of the sky, except for the 
area south of 5 = -85°, have been searched. In order to determine the 
galactic distribution from this sample of pulsars, one must take into 
account the selection effects of the surveys, which can be determined to 
reasonable accuracy. 

We use the method described by Taylor and Manchester (1977, hereafter 
referred to as TM) to calculate the distribution of pulsars as a function of period 
P, height above the plane Z, galactocentric radius R, and radio luminosity L. 

The data on 328 pulsars used here is from Manchester and Taylor (1980). The 
number of pulsars detected in the intervals P to P + dP, Z to Z + dZ, R to R + 
dR and L to L + dL will be: 

Nq (P,Z,R.L)dPdZdRdL = V(R,L)o(P,Z,R,L)dPdZdRdL 


( 2 ) 
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where V(R»L) 1s the effective volume of the GsIftKj/ between Z and Z + dZ 
and R and R + dR searched for pulsars with periods beitween P and P + dP 
and luminosities between L and L + dL. The true space density of pulsars 
In these same period and luminosity Intervals Is p(PtZ,R,L). As TM have 
argued, the distributions with respect to P and Z are not seriously altered 
by selection effects, and there seems to be little observational correlation 
between L and P. Therefore, V depends only on R and L. If the distributions 
in P, Z, R, and L are Independent, then the density may be expressed as: 

p(P,Z,R,L) = Pp(P)N(Z)D(R)pl(L). (3) 

The P and Z distributions, and N(Z), may be determined directly from 
the data without correcting for selection effects, while the R and L 
distributions will depend on V(R, L). 

a) Period Distribution 

There is evidence of a correlation between the P and Z distributions 
for pulsars in the second Molonglo survey (Manchester 1979, Taylor 1979). 
Pulsars with longer periods tend to have larger scale heights which would 
be expected if they are born with large space velocities from a population 
having a small scale height. While this correlation will affect the 
latitude distribution of the galactic y-ray flux from pulsars, it will not 
strongly affect the longitude profile, where we integrate over a latitude 
range which includes more than 2/3 of the pulsars. We will therefore 
treat Pp(P) and N(Z) as independent functions in computing the y-ray longitude 
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profile, but take this correlation into account in SV in computing the 
latitude profile. 

The period distribution for this sample of 328 pulsars is shown in 
Figure 1. The number falls off rather sharply at periods around 2s., 
which is well below the 3.9 s. long-period cutoff of all three surveys. 
The decrease in number of pulsars at short periods is also a real effect 
since the surveys were fully sensitive down to periods of .06s. The P 
distribution is therefore fairly narrow and peaks at periods between 0.5 
and 1,0 s. 


b. Pulsar Distances and the Interstellar Electron Density 


The Z artd R distributions depend strongly on the electron density 
distribution in the Galaxy, because pulsar distances are determined from 
their dispersion measures. If we assume an exponential dependence for the 
electron density. 


ng(z) = Oq exp(-|z|/hg), 


(4) 


then the distance to a pulsar with dispersion measure DM at galactic latitude 
b is given by 


d = 





(5) 


In the limit b -♦ 0°, equation (5) approaches d = DM/n^, so that the computed 
distance essentially varies inversely with n^. Since pulsars are a disk 
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population, local densities will vary approximately as n®. If one assumes 
that pulsars also have an exponential distribution, 

N(z) » Nq exp(-U!/hp), (6) 

then the pulsar scale height, hp will also depend strongly on n^ and weakly 
on the electron scale height, hg, as long as hg » hp. 

The best estimates of the average come from distance determinations to 
pulsars using measurements of the 21 cm. absorption of the pulsar signals by 
intervening hydrogen gas. From distance estimates for 32 pulsars and their 
dispersion measures, Weisberg et al. (1980) conclude that the mean inter- 
stellar electron density in the galactic plane on klloparsec scales away from the 
galactic center region is 0.02-0.03 cm-3. Although these is evidence for a 
longitude dependent variation in < 0 g> from this same data, we will assume that Is 
constant throughout the galactic plane. Within 1 kpc of the Sun, the dumpiness 
1n the electron density distribution becomes Important, and we take account 
of this by subtracting the contributions to the DM from nearby HI I regions 
in the line of sight according to the method of Prentice and ter Haar (1969). 

The scale height of the electrons is less well determined, but it must be 
greater than the pulsar scale height or there would be a cutoff in DMsih|b|, 
the Z component of dispersion measure, which is not observed. We will 
assume that hg = 1000 pc., which is the value adopted by TM, for the 
purpose of calculating distances. 

c. Z and R Distributions 

In light of the strong dependence of pulsar distances and densities on n^. 
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the estimates of which are somewhat uncertain, we have determined the Z and R 
distributions for several possible values of n^. Using the corrected dispersion 
measure da^' ind hg = 1000 pc., a least squares fit of equation (6) to the 
observed N(Z) distribution gives hp = 400 * 25 pc. for n^ = ,02 cm"®, hp = 325 
± 20 pc. for Og = .03 cm-® and hp = 230 * 15 pc. for n^ = .04 cm-®. These 
values are somewhat larger than the scale heights for n^ = .02 cm"® and .03 
cm”® derived by TM, probably due to the greater number of long period, high 
Z pulsars detected in the second Molonglo survey tyhich were not included in 
their sample. 

In calculating the distribution of pulsars in galactocentri'j radius, we 
have divided the sample into those with 0° < /• < ISO® (positive longitudes) and 
those with 180° sc / < 360° (negative longitudes). These two groups contain 
respectively 155 and 161 pulsars having measured fluxes at 400 MHz. We have 
used the iterative method described by TM to solve simultaneously for the R 
distribution and luminosity function for each group, assuming the 1 distribution 
of equation (6) and semicircular symmetry about the galactic center. The R 
distributions for n^ = .03 cm"®, in terms of the surface density of pulsars 
projected onto the galactic plane, are shown in Figure 2. The density scales 
were determined by integrating the derived luminosity functions to give the 
local surface densities. The most obvious feature at both ;’0s1tive and negative 
longitudes is a peak at around 5 kpc, where surface densities are five times 
greater than those at the Sun. The densities seem to fall off inside about 4 
kpc, but the statistics are very poor and no firm conclusions can be made about 
pulsar densities in the inner galaxy. Even though D(10) at negative longitudes 
is about twice as large as at positive longitudes, the maximum densities near 
5 kpc are very similar. The negative longitude distribution shows a 
secondary peak around 8 kpc and what appears to be a real deficit of 
pulsars between 6 and 7 kpc. These features are not present at positive 
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longitudes to any degree of significance. It is tempting to associate the 
secondary peak with the line of sight tangent to the Sagittarius arm which 
falls around 8 kpc. These distributions for = .03 cm“^ bear a striking 
resemblance to those of other Population I tracers. The 5 kpc peak in both 
ionized hydrogen (Lockman 1976) and CO emission, surveyed at positive 
longitudes only (Gordon and Burton 1976, Scoville and Solomon 1975), shows 
the same shoulder at 7.5 kpc and steep fall off at 4 kpc. as the pulsar 
distribution for 0 ^ X < 180°. The galactic y-ray emissivity, determined 
by unfolding the SAS-2 longitude data (Stacker 1977, Caraveo and Paul 1979), 
peaks around 5 kpc at both positive and negative longitudes but also has a 
prominent secondary peak around 8 kpc at negative longitudes, thus showing the 
same type of asymmetry as the pulsars. 

The R distribution changes drastically for other values of the mean 
electron density. The positive longitude distribution seems to be especially 
sensitive to the value of n^ assumed. For n^ = .02 cm"^ , the density increase 
toward the inner galaxy is much more gradual at both positive and negative 
longitudes. The 5 and 8 kpc peaks at negative longitudes shift inward and 
are less pronounced, while the structure at positive longitudes washes out 
completely, leaving no trace of a peak at 5 kpc. For n^ = .04 cm"^ , the 
density increase is steeper and the structure which is present at n^ = .03 
cm“^ again tends to wash out. At negative longitudes there is only one 
peak and it appears between 6 and 7 kpc. 

The greater sensitivity of the positive longitude distribution to 
changes in n^^, which are equivalent to changes in the distance scale, 
probably results from the way in which the observed pulsars are distributed. 
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At positive longitudes, the 5 kpc peak is mainly due to pulsars located 
between the Sun and the galactic center, where changes in distance result 
in large changes in R. Most of the structure in the negative longitude 
distribution is due to pulsars located along tangents to galactocentric 
circles, where changes in distance produce only small changes in R. Some, 
but not all, of this sensitivity to n^ at positive longitudes can be 
eliminated by removing the Arecibo pulsars from the sample. The Arecibo 
search was 10 times more sensitive than the others and confined to a small 
longitude range around t, = 40°. 

d. Numbers and Birthrates 

The total number of observable pulsars on each side of the Galaxy can 
be obtained by integrating over the R distributions: 

00 

Np = tt / D*(R)RdR. (7) 

0 

The numbers calculated in this way for the different values of n^ are listed 
in Table 1. One would expect that, since the Sun is not in any preferred 
position in the Galaxy, Ng and Ng would be equal. They turn out to be 
equal within the errors only for n„ = .03 cm"^ , with Nt < Np for n^ = .02 
cm"^ and Ng > Ng for n^ = .04 cm"^. This result would seem to argue in 
favor of a mean electron density of .03 cm“^ , and the argument is reinforced 
when one considers the shapes of the distributions for this value of n^ and 
their strong similarities to distributions of other Population I tracers in 
the Galaxy. 
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The total number of observable pulsars in the Galaxy for different 
n^ and the Sfdlactic birthrates which they imply are also listed in Table 1. 
The number Nq = (3.7 ^ 0.6) X 10® for n^ = .03 cm“’ agrees with the number 
Nq = (4.2 *1.6) X 10® obtained by K^nchester (1979), but is considerably 
greater than the result Nq = (1.3 -t 0.4) x 10® obtained by TM. If we 
assume that only 20 percent of pulsars are observable because of the 
beaming effect, and that their average lifetime is 10“^ year, which is 
probably an upper limit, then the galactic birthrate implied by this number 
is 1 every 5 years. Even considering the large uncertainties in this value, 
it would be inconsistent with the rate of occurrence of supernovae in the 
Galaxy if the rate were as small as 1 every 30 to 40 years (Milne 1979) or 
1 every 80 years (Caswell and Lerche 1979). 


\ 
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IV. LOCAl y-RAY PRODUCTION RATE 

Using the results of §II and 1 1 II, we can determine the local y-ray 
production rate for pulsars, 

qjsr^="100 MeV) = J L (P)pp(P)dP, (8) 

p o' 

where D'^^dO) is the surface density of pulsars at the Sun, hp is the scale 
height, [.^(P) is the y-ray luminosity formula given in equation (1) and 
Pp(P) is the period distribution in Figure 1, normalized so that 

/ Pp(P)dP = 1. 

0 

The value of q (> TOO MeV) can then be compared to the local production 

P W I 

rates which have been calculated for various diffuse processes. The 

ae 

luminosity L^(P) also depends on the surface magnetic field strength as B* , 
so that ppg^ (> 100 MeV) will be proportional to the average pulsar magnetic 
field strength. Since B cannot be measured directly, the only source of 
information on pulsar magnetic field strengths comes from measured slowdown 
rates and assumptions about neutron star structure. In view of this 
uncertainty, we will assume an average pulsar magnetic field strength of 
10^° gauss in making these calculations, but remember that the results 
could scale up or down with a change in this value. The production rates 
calculated according to equation (8) for the different values of n^ are shown in 
Table 2. The last column shows the fractional average production rate, <qpg^ '•.■q-j- 

~ (qpsr ^psr^^^%’ assuming a total local emissivity of q-j- = l.Sxin-^'^ cnr^ s"l . 
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The differential production spectrum above 100 MeV can be calculated 
from the pulsar spectra predicted by the model discussed in 5ll. The 
characteristics of these spectra were described in an earlier paper (Harding 
1981), where it was found that observed spectral data for both the Crab and 
Vela pulsars could be fit using the same parameters (except for the rotation 
period). If we assume that these parameters are valid for all pulsars, then 
their individual spectra will depend only on the period, and the differential 
production spectrum can be calculated using equation (8) with the intensity, 
(photons s"’ MeV’) substituted for the luminosity L (P) and 
<9pgr> in place of 9*5^. The resulting production spectrum for n^ = .03 cm"'’ 
is shown in Figure 3 along with production spectra which have been calculated 
by Stecker (1977) for the major diffuse processes. The calculated pulsar 
spectrum has been extrapolated below 100 MeV, although the model may not be 
accurate at the low energies due to the undetermined synchrotron contribution 
from secondary particles. The extrapolation thus represents a lower limit on 
the pulsar contribution below 100 MeV. Above 100 MeV and below 1 GeV, the 
pulsar spectrum is almost as steep as the bremsstrahlung and Compton spectra, 
and becomes steeper above 1 GeV due to the pair production and curvature 
radiation cutoffs. The pulsar production rate exceeds the "standard" 
production rate for electron bremsstrahlung calculated by Stecker up to 
about 2 GeV, and thus provides an additional source of emission which can 
steepen the total spectrum (Harding and Stecker 1980). 

This steeper total spectrum is a better fit to the observed spectrum 
of the galactic emission than is the total spectrum from only diffuse 
processes. Figure 4 shows the total production spectrum of Figure 3 and 
the pulsar spectrum, which have been adjusted arbitrarily in magnitude to 
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I 

fit the data points. There is no evidence in the data for a prominent 
peak at 70 MeV from tt° decay, though the points cannot be fit by a straight 
power law either. It is interesting to note that the pulsar spectrum alone 
fits the data quite well, so that there are no spectral constraints in this 
model on the size of the pulsar contribution. In fact, pulsars could be a 
major source of the y-ray emission below 100 MeV. 
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V. LONGITUDE AND LATITUDE PROFILES 

The distribution of pulsar Y-ray flux in the Galaxy as a function of 
longitude and latitude may be determined from the local production rate 
and the galactic distributions, determined in §III. The flux from a given 
area of the sky dxdsinb will be 


d5(x,b)=5£sjd J Pn(s,x,b) ds dsinb dx, (9) 

4tT 0 

where s is the distance along the line of sight and pp(s,X,b) is the pulsar 
density. 


The flux at each longitude, integrated over a certain latitude range 
b ^ b„^„, gives a longitude profile: 

liiaA 


d§(jk X) = ‘^psr 
dx 4'H’ 


bmax „ 

i* .]■• 

“^max ° 


Pp(s, ± X, b) ds dsinb. 


The pulsar density is 


( 10 ) 


Pr(s. ± X. b) = » , ^ .^1 NU(s,b)}, (11) 

D=^(10) 

where N(Z) is the exponential Z distribution given in equation (6) and D'*'(R) 
is the R distribution. The coordinate transformations between (R, Z) and 
(s, X, b) are 


Z = s sinb, 
R 

r = s cosb. 


(R| + r^ 


- 2Rg r cosX)‘ 


( 12 ) 
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where r is the distance s projected onto the plane and Rg = 10 kpc is the 
distance of the Sun from the galactic center. If these expressions are used 
in equation (11) to determine the density at each point along the integration 
path, then equation (10) becomes: 

d$(±X) = W h J D^rR(r, /)} Cl - expC*" Vx/h )} ~ (13) 

aT 0 D*( 1 0) 

We have numerically integrated this expression to obtain the y-ray 
flux in 10° longitude bins over a latitude range (b| s: = 10°, using 

a step size in r of 200 pc. The resulting profiles for the three values 
of electron density, n^ = .02, .03, and .04 cm"^ , are shown in Figure 5. 

The flux level increases as the mean electron density, and therefore the 
pulsar density, increases. The n^ = .03 and .04 profiles show sharp 
increases in flux between X = 300° and X = 60°, with peaks in emission around 
X = 30° and X = 330°. They are roughly symmetric around X = 0°, although 
the flux decreases toward the anticenter more slowly at negative longitudes. 
The n^j = .02 profile is very asymmetric, with a much greater flux at negative 
longitudes. The large asymmetry in the R distribution for n^ = .02 was 
discussed in §III and it indicates that this value of n^ is probably too low. 

These calculated profiles are plotted with the SAS-2 lo'^gitude data 
in Figure 6. The contribution from identified 'point sources has not been 
subtracted, but the strongest of these sources, including the Crab and 
Vela pulsars, are labeled. The calculated profiles are similar in shape 
to the observed profile , but they have a larger center to anticenter 
flux ratio. For n^ = .03 cm"^ , the pulsar contribution to the total flux 
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is 15 percent in the interval 300° < < 60° which is about twice as large 

as the 7 percent contribution in the interval 60° < * < 300°. For n^ = ,04 
cm"^, the center contribution goes up to nearly 25 percent, while the anti- 
center contribution is around 10 percent. As mentioned in the last section, 
the pulsar y-ray production rate scales with average magnetic field strength. 
Therefore, for an average surface field of say, 2 x 10^® gauss, the fluxes in 
Figure 6 would all ^e twice as large. 


The latitude distribution of pulsar y-ray flux can be calculated for a 
particular longitude or longitude range. Using equations (9), (11), and 
(12), we have 


/.»rL,(P)Pp(P) ^ip^«pC-s|sinb|/hp(P)ldPds. 


(14) 


Because of the dependence of scale height, h , on period as discussed in 

r 

5lII, the production rate is also a function of P and must appear under the 
integral along with the distribution functions. A least squares fit to the 
plot of average Z distance versus period for the pulsars in the Molonglo 
survey (Figure 2 of Taylor 1979) gives: 

0.3E4 

hp(P) = 396 P pc. (15) 


The Z distances were calculated assuming that n^ = .03 cm"^ and hg = 1000 pc. 
Using this relation in equation (14) produces the latitude distribution for 
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X * 0° shown in Figure 7. The profile is quite narrow due to the fact 
that the majority of y-rays are produced by short period pulsars which 
have a relatively small scale height. At i » 0°, most of these pulsars 
are at a distance of around 5 kpc, producing a distribution with a 2° 
half-width. Because of the dependence of scale height on electron 
density, this width would be larger for n^ = .02 cm"^ and smaller for n^ 
= .04 cm"^. Since the pulsar latitude profile is so narrow, the size of 
the contribution to the total flux is not restricted by the small width 
of the observed latitude profile. 





VI. DISCUSSION 
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We have presented a number of new results in this paper both on the 
distribution of pulsars in the Galaxy and on the pulsar contribution to diffuse 
galactic y-rays. The pulsar R distribution 1s found to be significantly 
different at positive and negative longitudes. This difference may result 
from the spiral structure of the Galaxy if pulsars are located preferentially 
along spiral arms. The sensitive dependence of the distributions on mean 
electron density in the galactic plane provides a strong argument in favor 
of Oq » .03 cm*5. This value is needed to give equal numbers of pulsars on 
both sides of the Galaxy and distributions which resemble those of other 
Population I objects. 

The y-ray flux contribution from pulsars in this model depends both on 
the number of pulsars in the Galaxy s^d on their average magnetic field 
strength. The two largest sources of uncertainty in the number of pulsars 
(as determined from the local surface density) are the mean electron density 
in the galactic plane and the minimum pulsar luminosity. The range in electron 
density we have considered, .02 cm"^ d n^ s .04 cm-®, represents about a factor 
of 3 uncertainty in the pulsar flux contribution. However, the arguement in § III 
(d) seems to rule out densities as low as .02 cm"®, so the range of possible values 
of Oq is pV-obably smaller. There is evidence that the electron density distri- 
bution is much more complex than the simple exponential model used here. The 
densities seem to be significantly higher in the inner Galaxy (Weisberg et aU 1980), 
suggesting an R dependence, and higher south of the galactic plane, suggesting a * 
multicomponent Z distribution (Harding 1980). A model for the electron density 
incorporating these features could change the derived pulsar densities (Manchester 
1980). 

The local pulsar density derived for each value of n^ depends critically on 
the lower limit chosen on the integral over pulsar radio luminosities. Due to 
the limited sensitivity of the searches, there are fewer low luminosity 


pulsars in the sample and thus, some uncertainty in where the derived power 
law luminosity function turns over. However, a turnover is expected somewhere 
around the present sensitivity limit because the existence of many more 
pulsars in the Galaxy would only widen the gap between the pulsar birthrate 
and the supernova rate (cf, § III d). 

Pulsar magnetic field strenvjths can be indirectly inferred from their 
slowdown rates. Although the rates of slowdown can be measured fairly 
accurately for many pulsars, a derivation of the surface magnetic fields 
requires a knowledge of neutron star radii and moments of inertia, which 
can only be determined by knowing the correct equation of state. The 
present theories of neutron star structure include a wide range of possible 
equations of state which predict a range of stable masses and radii (Baym 
and Pethick 1979), The range in these parameters gives a factor of 2 to 3 in 
average magnetic field strength, which gives the same uncertainty in the 
pulsar flux contribution. The range of average field strengths for a neutron 
star mass of 1,4 Mg is roughly (0.7 - 1.6) x lO’* gauss, with the lower 
fields more likely for stiffer equations of state. The local y-ray production 
rate from pulsars for n^ = .03 cm"^ and an average field strength of 1 x 10^® gauss 
is 15 percent of the total production rate (cf. Table 2). This fraction 
could be anywhere from 10 to 25 percent from the above field strength range. 

There are several characteristics of the galactic pulsar y-ray flux 
which emerge regardless of the uncertainties. The production spectrum is 
as steep as the bremsstrahlung and Compton spectra above 100 MeV, although 
it is not a straight power law. The fact that the shape of the pulsar 
spectrum fits the observed galactic plane y-ray spectrum quite well allows 
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the pulsar flux to be any fraction of the total emission. The shape of a 
pure ’T® decay spectrum does not fit the COS-B spectral data points (Paul 
et al. 1978), placing an empirical upper limit to the fraction which 
decay can contribute to the total emission. Several authors (Cesarsky 
et al. 1978, Strong et al. 1978, Hartman et al. 1979) have concluded that 
large bremsstrahlung contributions are needed to give total spectra which 
fit the shape of the observed y-ray spectrum. The size of these contri- 
butions requires a flux of cosmic ray electrons which 1s from 2 to 5 tliries 
higher than what Is expected from the observed electron flux at the Earth, 
after correcting for solar modulation. In addition, the ratio of bremsstrahlung 
flux to decay flux needed requires much higher electron to proton ratios than 
are observed near the Sun, Pulsars provide an additional steep spectrum source 
of emission which has not previously been included In analyses of the galactic 
Y-ray spectrum. With the pulsar contribution Included, a good fit to the 
spectral data can be obtained with a lower bremsstrahlung contribution and 
thus, lower electron fluxes and electron to proton ratios. 

We have also found that the pulsar flux percentage contribution is twice 
as large toward the galactic center as it Is toward the anticenter. The 
diffuse flux component, which is the difference between the total flux and 
the point source flux, will then haVf a smaller center to anticenter ratio 
than the total flux. The cosmic ray density variation in the Galaxy needed 
to explain the total flux Increase toward the galactic center could therefore 
be smaller. The latitude distribution of the pulsar flux is quite narrow 
and, in fact, has a smaller width than the observed latitude distribution 
of T-ray flux toward the galactic center (which may be resolution limited). 

We have obtained this result by using recent direct evidence that short 
period pulsars, predicted to be the most luminous y-ray emitters, have a 


\ 
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smaller scale height, comparable to their Population I progenitors. Studies 
of the diffuse galactic /-ray spectrum at high latitudes (Fichtel et al. 1978; 
Lebrun and Paul 1979) also show the significant low energy excess seen at low 
latitudes, requiring a steep spectral contribution greater than the “standard" 
bremsstrahlung component. From Figure 7, the estimated pulsar flux is about 
10 percent of the galactic flux at b ~ 15* , given by Fichtel et al. to 

be (4.0 ± 1.0) X 10"® cm"* s”’ sr-’ , Even at these higher latitudes, then, 
pulsars could supply some of the steep spectral contribution that is needed. 
This high latitude contribution would come primarily from older pulsars, 
individually weak as y-ray sources, and therefore consistent with the absence 
of strong high latitude point sources. 

In estimating the pulsar contribution, we have counted only observable 
radio pulsars as sources of y-ray emission. It is possible that there exists 
a population of pulsars emitting pulsed y-rays but no pulsed radio emission. 
In the y-ray pulsar model considered in this paper, the shortest period 
pulsars have the highest luminosities. If there are long period pulsars 
which turn off in the radio before turning off at y-ray energies, they 
would not be luminous enough to contribute much to the galactic emission. 

If, however, there are short period pulsars which have turned on at y-ray 
energies prior to turning on in the radio, then these objects could be 
contributing as much as the radio luminous pulsars to the galactic emission. 
Michel (1978) predicts that incoherent radiation makes up an increasingly 
greater fraction of the total available energy as pulsar period decreases. 
Since observable radio emission requires a coherent process, pulsars with 
extremely short periods could be under luminous in the radio but very 
luminous in y-rays. 


In summary, the present estimate of a 15 percent to 20 percent pulsa" 
contribution to the total y-ray emission in the Galaxy represents only ’the 
contribution from conventional radio pulsars. It indicates that the total 
contribution from all point sources of y-ray emission could be substantial 
and should not be overlooked in considering models for diffuse processes and 
cosmic rays. 

The author would like to thank F. W. Stecker and D. J. Thompson for 
extremely useful discussions and suggestions, and J. H. Taylor and 
P. R. Backus for supplying up-to-date pulsar data. 
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Figure 1 
Figure 2 

Figure 3 


Figure 4 
Figure 5 
Figure 6 
Figure 7 


FIGURE CAPTIONS 


--Period distribution of the observed pulsars. 

—Surface density of pulsars versus galactocentric radius at 
positive and negative longitudes for a mean electron density, 
n^ = .03 cm"®. 

—Differential production spectrum for pulsars with an average 
magnetic field strength of lo'® gauss, for a mean electron 
density n^ = .03 cm-®. Production spectra for diffuse processes 
are from Stecker (1977). The high energy part of the rP 
spectrum is from Stecker (1979). 

—Total spectrum and pulsar spectrum from Figure 3 fit to the 
data points from SAS2 and COSB for 355° < JL < 15°. 

—Pulsar y-ray flux longitude profiles for different values of 
the mean electron density in cm"-®. 

—Longitude profiles from Figure 5 plotted with the SAS2 
longitude flux from Hartman et al. 1979. 

--Latitude profile of pulsar flux at A = 0° for n^ - .03 cm"®. 






CM 

I 

u 

Q. 

(O 

(T 

< 

{/) 

-J 

ZD 

Q. 


2000 

1500 

1000 

500 

2000 

1500 

1000 

500 






Q(E^)( PHOTONS cm s MeV 




PHOTONS cm s rad MeV 


10 


-5 


10 


rS 


10 


-7 


lO’® h 


TOTAL 



f SAS 2 (HARTMAN et al 1979) V 
COS B ( PAUL et al 1978 ) 

^ L 

10 ^ 10 ^ 


10 


Ey(MeV) 


Figure 4 







VELA 

(PSR0833-45) 



(QVa-S- WO)/(A®WOOI<3)SNOiOHd 
2 


\ 


Fif.ure 


01 X) 




o 


(js-s-,tuo)/SNOiOHd 


i 


Figuy'e 7 
\ 


o 

Csj 


